arXiv:l504.05219vl [math.ST] 20 Apr 2015 


BIVARIATE NATURAL EXPONENTIAL FAMILIES WITH 
QUADRATIC DIAGONAL OF THE VARIANCE FUNCTION 


JOANNA MATYSIAK 


Abstract. We characterize bivariate natural exponential families hav¬ 
ing the diagonal of the variance function of the form 


diagL(mi, m 2 ) = (Am? -|- ami + hm 2 + e, Am| -|- cmi -|- dm 2 -I- /) , 


with A < 0 and a, ...,/€ R. The solution of the problem relies on find¬ 
ing the conditions under which a specihc parametric family of functions 
consists of Laplace transforms of some probability measures. 


1. Introduction 


Let us recall some basic concepts concerning natural exponential families 
(NEFs for short). For a positive measure n on M” we define its Laplace 
transform by 



Let 0(/t) = lnt{0 € M"" : L^{9) < 00 } and let kfj, = logL^ denote the 
cumulant function of Set M. as the set of such measures fj. that 0(/i) 7 ^ 0 
and fi is not concentrated on an affine hyperplane of M"'. For /r € M, the 
family of probabilities 

F(/i) = {P{0, ^i){dx) = exp {< 9,x > —k^{9)) fj,{dx) : 9 € 0(/i)} 

is called the natural exponential family generated by ji, see [15]. The image 
Mp of 0(/t) by diffeomorphism k'^ is called the domain of means of fi. The 
variance funetion (VF) of the NEF is defined by yF(m) = A:"(T^(m)), 
where m G Mp and = {k'^~^. In what follows we will be frequently 
omitting the subscript F in the notation Vp for the variance function. 

The principal significance of the mapping m 1 —>■ I/p(m) is that, together 
with its domain Mp, it characterizes NEF uniquely. It makes it possible to 
describe NEFs assuming a concrete form of the variance fnnction. There has 
been a lot of interest in such questions. For example, Letac m characterizes 
all NEFs with VF(m) = Bm + C, where R is a linear operator mapping 
M"" into (the space of n x n symmetric real matrices), and C G Sn- 
Casalis |4], [5] gives a generalization of this result by considering I/p(m) = 
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am (g) m + i?m + C. Her result was further extended by Hassairi and Zarai 
[ 8 ] to NEFs with cubic variance functions. 

For another standpoint, we recall some papers that provide characteri¬ 
zations based on a fragmentary knowledge of the variance function. Kok- 
endji and Seshadri [TO] start with det V (m) = const and thus identify the 
Gaussian law in M”. In [14] Letac and Wesolowski characterized NEFs with 
V (m) = p~^in (g) m — where M^, is a symmetric matrix associated 

with a quadratic form m —)■ (;/!)(m) an unknown real function, and p is a 
number. 

One of the most important papers that base on a partial knowledge of Vp 
is [T], in which the authors considered the diagonal family of NEFs in M” 
such that 

diagl/(m) = (/i(mi),...,/„(/«„)) 

(diagF stands for the diagonal elements of the variance matrix V, and /* is 
an arbitrary function of i-th coordinate of m). They gave the full charac¬ 
terization of the family (showing also that /,, i = 1,2,... ,n must be some 
polynomials of degree at most 2 ). 

A natural question is the following: let us fix n classes of functions Ti 
of the variables m = (mi, m 2 ,..., m„) for f = 1,... , n. For instance Ti 
could be the functions of the form /i(m) = Amf + B{m) + c, where H is a 
linear form and A, c € M. What happens if we stay focused on the diagonal 
(/i(m),..., /n(m)) of the VF such that fi € Ti! We are going to study this 
question for n = 2. For a partial result in this direction, see |6|. 

The paper is organized as follows. Its main results. Theorem 12.21 and 
Theorem Ea are formulated in Section [ 2 j Sections [3] and [H provide proofs 
of the theorems, along with some auxiliary facts. The straightforward, but 
lengthy proof of one of the facts. Proposition 14.21 is presented in Appendix 
A. 


2. Main Results 

Our task in this paper is to characterize NEFs with the VF of the form 

( 2 . 1 ) diagV(mi, m 2 ) = {^Am\ + ami + bm 2 + e, Am^ + cmi -|- dm 2 + /) , 

with A < 0 and a, ...,/€ M. To avoid considering the VFs of the diagonal 
NEFs analyzed by Bar-Lev et al. in [1], we need to exclude the case b = 
c = 0. Due to the symmetry between b and c, in the sequel we shall assume 
without loss of generality that b ^ 0. 

We decided to stick here to A < 0 only and to place our findings on the 
case A > 0 in a separate future paper. Our decision is motivated by a rather 
surprising lack of symmetry between the two cases. It turns out that if A > 0 
then ( 12 .Ih implies a much wider and less homogeneous class of parametric 
bivariate functions that are candidates for the Laplace transforms of the 
corresponding NEFs, than for A < 0. 
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Remark 2 . 1 . An interesting feature of the problem is that it can be formu¬ 
lated equivalently using regression properties. Assume that X = (Xi, A 2 ) 
and Y = {Yi,Y 2 ) are independent identically distributed random vectors 
such that their Laplace transform exists in a set with a nonempty interior. 
The regression formulation of the problem we solve in this paper is to iden¬ 
tify the distributions of X if 

E((Xi-yi)2-2AXiyi|X + Y) = aiXi + Yi) + b{X 2 + Y 2 ) + 2e, 

E((X2 - Ya)^- 2 AX 2 Y 2 IX + Y) = c{Xi+Yi) + diX2 + Y2) + 2f. 

For other works in a similar vein, see e.g. m, m, i, m, m- 

Now we state the main theorem, that specihes the form of a Laplace 
transform corresponding to the measure generating NEF with the VF given 
in ra . 

Theorem 2.2. Let ^ be a probability measure (not concentrated in a point) 
generating NEF with VF given by m- Then the Laplace transform of p, 
is of the form 
( 2 . 2 ) 

exp (^Vi +, ( 01 , 02 ) G 0 (h), 

where {Ai,..., is the set of distinct real roots of 

(2.3) A^ — 2aA^ -I- (2Ae + — db) A^ — {2Aae — adb + ch^) \+ 

- edbA + fb^A = 0 

(the coefficients in (j2.3jl come from ()2.ip ). Furthermore, Ur > 2 and either 

• oti > 0, i = 1,..., Ur, —1/A € N and 0(/r) = or 

• cii < 0, i = 1 ,..., Ur, ~ “1/^ ^ 2N and Q{p) = 

Remark 2.3. It is clear that if one of the bullet item conditions above holds 
then ()2.2p is the Laplace transform of a probability measure. 

It might be worthy to emphasize that the only assumption we make about 
0 (/r) is that it contains a neighborhood of the origin (this follows from the 
fact that p belongs to NEF). The concrete form of Q{p) {&{p) = E^) is not 
assumed but it is implied by the assumptions of Theorem 12.21 

The proof of Theorem 12.21 contains two parts. The first part explains why 
the distribution sought in (12.Ih . must have the Laplace transform of the 
form (12.2p . We present this partial result in Proposition 14.21 The second 
part consists of some arguments concerning the conditions on A, a'^s and 
&{p)- The arguments are gathered in a stronger result (see Theorem 12.41 
below), which, by a special choice of a matrix A in the case of Ur > 3, 
implies the conditions on and A (the case = 2 will be considered 
separately in the proof of Theorem 12.2h . 

The abovementioned stronger result is the following 
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Theorem 2.4. Let r > 0 , A = (Ai,A 2 ,A 3 )'^ G and suppose 0 C 
contains a neighborhood of the origin. Define 

(2.4) L{9) = ^ao + aj exp (At, 9)^ , 0 G 0. 

Assume that equation 

(2.5) a'^A = 0 

has no solutions in a = ( 01 , 02 , 03 ) G such that 3i^j UiOj < 0. Then L is 
Laplace transform of a probability measure if and only if either 

(a) ai > 0, i = 0,... ,3, = 1 and r G N, or 

(b) Oi < 0, i = 0,..., 3, = “1 and r G 2N. 

Notice that conditions (a) and (b) are sufficient for L to be the Laplace 
transform of a probability measure. 

Remark 2.5. Bar-Lev et al. [Proposition 3.1(a), [T]] consider a similar class 
of Laplace transforms. If n G N, define T to be the family of non-empty 
subsets of { 1 , 2 ,..., n}, and for z = (zi,..., Zn) put 

z'^ = Zj, for T eT. 

j&T 

Proposition 3.1 from [1] describes for which a the analytic function of z 

(2.6) (1 + AT G N, 

Ter 

has non-negative coefficients in the Taylor expansion in the neighborhood of 
the origin. For Zj = e^fi j = 1,..., n, this is equivalent to deciding whether 
(12. 6 p is the Laplace transform of a positive measure. 

Observe that our approach differs from the one taken in [T]. This is because 
in our problem (|2.4I) . the crucial step is to justify the series expansion of 
L, while in (|2.6p analyticity is assumed. The A matrix and the support 
structure condition (j2.5ll allow as to identify the masses of atoms among the 
series coefficients and to infer about the signs of a’s. As a consequence of 
the need to justify the series expansion, we restrict ourselves to Another 
difference lies in the fact that the support of the measure corresponding to 
(j2.6ll belongs to a lattice in The support of the measure corresponding 
to ( 12 .4p depends on the choice of A. 

Remark 2.6. Again, it might be worthy to emphasize that our assumptions 
on 0 in Theorem 12.41 are weak. They are dictated by the fact that the 
measure with Laplace transform (|2.4p belongs to a NEF. 
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3. Proof of Theorem 12.41 


Proof of Theorem 12-41 Let us consider first L on 'D: 

= i 0 G Qfj exp (Aj, 0) > 0 

I i=0 

We shall be using the following relation between two vectors x = (xi, X 2 , X 3 ) 
and y = (yi, 2 / 2 , 2 / 3 ) in 



X xi < yi V (xi = ?/i A X2 < ^2) V (xi = yi A X2 = y2 A X3 < y^). 

Observe that (12.51) implies the existence of permutation a of {0,1, 2,3} such 
that 

(3.1) -^ITO ^ P A(J 2 P Act 3 , 


where Aq = ( 0 , 0 , 0 ). We will split the reasonings into four cases. 

(1) Assume first that auQ > 0. We will show that cio-. > 0, f = 1,2,3 
and r € N. Note that (j3.ip yields 


(3.2) 


An A„, - A^, ^ A^, - A„r, p - A„ 


This in turn allows us to choose 9* £ P (negative with arbitrarily 
large modulus) to ensure that 
3 


^^exp(A^, -A^^,e*) 


< 1 . 


Define H{6) = exp (—rAg-g, 0) ^ ■ {H is the Laplace transform 

of a probability measure if and only if L is the Laplace transform of a 
probability measure.). Function H has the following series expansion 
in the neighborhood of the origin: 


H{6) = cY^ 

j=0 


r{r-!)■■■ {r - j+ 1) 


— exp (A,^, - A^g ,9* + 9) 


.1=1 “'^0 


where C = (^ exp (A^, - A^g ,9*)) ^ 

From (|3.2p we see that C'rL^exp(ACTi — Ao-g,0*) is the only coef¬ 
ficient at exp (Ao-^ — Ao-g, 9) in the expansion of H, hence > 0. 
Since (j2.5ll assures that there is no a G N, such that A^-^ — Ag-g = 
o (A 0-3 — Ao-g), we get that Crexp (Aq-^ — Ag-g, 9*} is the only coef¬ 
ficient at exp (Ao -2 — Aq-q, 9) in the expansion of ff, so — 0- Anal¬ 
ogously, from (j2.5[) there are no oi, 02 G N such that A^-g — A^-q = 
ai (Afjj - A^g) -F 02 (A<j 2 - Afjg), hence Cr^ exp {A^^ - A^^g, 0 *) is 
the only coefficient at exp (Ag-g — A^-g, 0) in the expansion of H, and 

cr^jg P 0 . 

Now suppose that r ^ N. Consequently, there exists the small¬ 
est k such that r(r — 1) • • • (r — /c -|- 1) < 0. Since coefficient at 
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exp [k (Ao-i — Ao-q) , 9) is positive as a mass of an atom in k (Ao-j — Ao-p) 
there exist ni, n 2 , ^,3 € N such that k (Ao-i — Ao-p) = ni (Ao-i — Ao-p) + 
^2 (A(J 2 — Ao-o)+n 3 (A 0-3 — Ao-p) which contradicts (j2.5[l : hence r E N. 
(2) Now assume that Ouq < 0 and > 0. We will show that in fact 
“(TO = Oj CKci > 0, z = 1,2 and r E N. 

Note that (l.l-ip yields 

(3.3) .^(70 ^ ^cr\ ^crz ^ ■^(72 -^(73 ^ Ag. 


Therefore we can choose 6* ^T) satisfying 


E 


j=0 


^exp(A,^. 

LX(j^ 


A<,3,r) 


< 1 


and define G{9) = exp {—rA^, 9) ^ ■ [G is the Laplace transform 

of a probability measure if and only if L is the Laplace transform of 
a probability measure.) We can write its series expansion: 


G{9) = gY, 

j=0 


r{r — 1 ) • • • (r — j + 1 ) 


^ — exp (A^, - A<^ 3 , 6 » + 9*) 


i=0 


where G = exp (A^, - A ^3 ,0*)) . 

From (13.3p we see that Gra~^acr 2 exp (Ao-j — Au^,9*) is the only co¬ 
efficient at exp (Ao-2 — A0-3, 9 ) in the expansion of G , hence 00-2 > 0. 
Using the same reasoning as in the preceding case, we conclude 
from (j2.5p that C'ra“ 3 ^ao-i exp (Ao-^ — Aa-^,9*) is the only coefficient 
at exp {Aai — Aa- 3 ,0) in the expansion of G , hence Oo-i > 0, f = 0,1, 

so Oo-p = 0. 

In order to conclude that r E N it suffices to repeat the previous 
reasoning. 

(3) Now assume that < 0 and cto-p < 0 (or < 0 and cto-p < 0). 
We will show that this assumption leads to contradiction. 

Since a^i + 00-2 = 1 “ Q^o-o “ conclude that a^i > 5 or 

0(72 > Without loss of generality we let Oo-j > 

Again (|3.ip implies 


Aao Act]^ ^ Aq A(J2 Aq-]^ ^ .^(73 A(j ^. 

Set T{9) = exp (—rAo-i, 9) L{9) and expand it in a neighborhood 
of the origin: 


T(9) = £ 


r(r — 1 ) • • • (r — /c -|- 1 ) 


fe=o 


kla^_^ 


E 


a. 


r, exp {Aai - A^i, 9) 


From (|2.5p . rarjia^~^ is the only coefficient at exp {Aa■^ — A^j, 9) in 
the expansion of T, hence Oo-. > 0 for z = 0, 2,3, so a^o = ®t 73 = 0. 
It contradicts the assumption. 
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(4) Now assume that = CKo-g = 0. A simplification of the analysis 
given in cases (1) and (2) leads to > 0 and > 0. Once we 
are done with it, the fact that r € N follows from the same steps as 
in the first case. 

Summarizing, we have shown that L (defined on D) is a Laplace transform 
of a probability measure if and only if (a) holds. 

Assume now that r € 2N. If so, we can consider L 

L(6»)= Q!o + ^aiexp(Ai,6') = -oq - ^ exp (A^, 6») , 


i=l 


i=l 


on 


V = ieG 


: ^ajexp(Aj,6') < o| = 

i=0 J 


9 G 


: ^-ajexp (Aj,6») > o| . 

i=0 ) 


Dehning dj = —a*, i = 0,... ,3, from (a) we can conclude that d, > 0, 
hence a, < 0, i = 0,..., 3. Thus we arrive at (b). □ 

4. Proof of Theorem 12.21 

4.1. Prom the form of the diagonal (j2.ip to the characteristic equa¬ 
tion (12.3p . Here we shall justify the transition from the diagonal (12.ip to 
the characteristic equation (j2.3p . 

Let A: be a cumulant function of a measure /i belonging to NEF satisfying 
(12.11) (that is k{ 9 i, 62 ) = logL( 0 i, 02 )j (^ 1 ,^ 2 ) S ©(z^))- Then condition (j2.ip 
can be written equivalently as 

d'^k . f dkV dk 8k 


(4.1) 


(4.2) 


del 

d‘^k 


= A 


del 


dk \ 


dk 


dk 


_ O — ^ - ^ ^TUi-^ f 


dej 


(901 502 


Define R = e . Then (14.ip and (14.2p become 


(4.3) 


(4.4) 


d'^R 

d'^R 

Wl 


dR dR 

SR , MR 


Since we assume that 6 7 Z Q, as a consequence of (14.311 and (14.411 we get 
d^R ^ d^R . , 2 /n . 9R 

+ — edbA + fb^A = 0, 


(4.5) 
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with the characteristic equation (j2.3p . 

Remark 4.1. If instead of making the assumption on b, one assumes that 
c 7 ^ 0 , then, due to symmetry between b and c, one obtains 

+ A^f'^ — acfA + ec^A = 0, 

which is analogous to (j4.5p . with the characteristic polynomial 
v^—2dv^+{2Af — ac + (f) (bc^ — acd + 2df A) f'^—acf A+ec^ A = 0. 

The roots of this equation in and the roots of (|2.3|) are connected via 

(4.7) = aA + — eA, 

(4.8) = c\ + du — fA. 

Careful analysis of the roots of (j2.3p leads, via various solutions of (j4.5p . 
to the following 

Proposition 4.2. Let fi be a probability measure (not concentrated in a 
point) generating NEF with VF given by (12.ip . Let {Ai,..., A„^} be the set 
of distinct real roots of (12.3p . Then Ur >2 and the Laplace transform of // 
is 

(4.9) L{9i, 92) = ai exp ^Ai^i + ^^——92 
with ( 01 , 6 * 2 ) G &{r)- 

The proof will be presented in Appendix lA. 11 

Remark 4.3. Analyzing the form of the Laplace transform L given in Propo¬ 
sition 4.2, one can see that the measure corresponding to the case A = —1 
consists of point masses which are concentrated on a parabola (see (14.71) 1. 

Now, we are in a position to prove our main result. 

Proof of Theorem \2.2[ Let (X,Y) be a random vector with Laplace trans¬ 
form (1221). If 

(X, Y) = {X - Ai/A, bY + aX + eA- A?/A) 
then the Laplace transform of (X,y) is 

( Ur 

ai + '^ai exp [(A^ - Ai) 0 i {\} - Ai) 02 ] 
i=2 

where (0i,02) G ©(/I) and ft = C{X,Y). We shall be working with the 
Laplace transform ()4.in|) rather than with ()2.2h because the form of (I4.10p 
allows us to use Theorem Ea 
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Now we will split our reasonings with respect to the number of distinct 
roots of (12.3p . 

First, if = 3 or 4, we will indicate how to specify A in Theorem 12.41 
to make it answer the questions about the coefficients from (I4.ini) (as a 
consequence also for (j2.2[) l. 

• For = 4 we choose 




A2 — Ai 

4-4 

0 ■ 

(4.11) 

A = 

A3 — Ai 

4-4 

0 



_ A4 — Ai 

4-4 

0 


Obviously this particular A plugged into (12.4h yields (j4.10p (after 
changing aoj «!) “2, ct3 into ai, q; 2,03,04, respectively). What we 
need to show is that snch A satisfies (j2.5p . To this end we first 
assume without loss of the generality that A2 < A3 < A4. Eqnation 
(j2.5p can be rewritten equivalently as the system of three eqnations 

fli(Ai — Ai) + a2(Af — A3) = 0, z = 2,3,4. 

From the first two we get 


fli + ci2{Xi + Ai) — 0, i — 2,3. 

Subtracting the first {i = 1) from the second one (i = 2) yields 


a2(A2 — Ai) — 0, 


hence 02 = 0 and ai = 0. Therefore (14.lip satisfies (j2.5p . 
• For Ur = 3 we take 


(4.12) 


A = 


A2 — Ai A2 — A^ 0 
A3-A1 Xl-Xj 0 
0 0 0 


Analogonsly to the preceding case, after plugging (I4.12p into (12.41) 
we obtain (|4.10p and the same analysis shows that (|4.12l) satisfies 

(HAp. 

For Hr = 2, (|4.10l) becomes 

L(6(i, 02) = (ai + 02 exp((A2 - Ai)0i + (A 2 - A?)02)) ^ • 

In order to prove the necessary condition, it is enough to analyze one¬ 
dimensional Laplace transform I given by 

l{6) = L{9/ (A2 - Ai), 0) = (oi + 02 exp(0))“^ 


in a neighborhood of 0. Our aim is to show that the necessary conditions 
for Z to be a Laplace transform of a probability measure on M, are either 

(1) oi, 02 > 0, oi + 02 = 1 and —^ € N, or 

(2) oi, 02 < 0, 01+02 = —1 and —^ G 2N 
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(sufficiency of the conditions is clear). 

Let us consider I onVi = {9 : ai+a 2 exp(0) > 0}, what implies ai+a 2 = 
1 (recall that /(O) = 1). Therefore, without loss of generality one can assume 
that > 1/2 (so a 2 < 1/2). What we want to prove is that ai > 1/2 
implies 02 > 0 . If ai = 1/2, then 02 = 1 — cti = 1/2, so let us assume 
that ai > 1/2. In such case |a2/ai| = |(1 —oi)/«i| < Ij hence, in a 
neighborhood of the origin, I can be written as 
(4.13) 


m = a, 2 ^- - - 

k=0 



k 

exp{k9). 


The only coefficient at exp(0) in (|4.13l) is a 20 L^ hence a 20 i-^ ^ > 0 

and so a 2 > 0. 

Now we focus on the part of conclusion dealing with the exponent r = 
— IjA. We provide here a reasoning analogous to the one given for r in the 
proof of Theorem 12.41 Suppose that r ^ N. Then there exists the smallest 
integer k such that r(r — l)...(r — A; + 1) < 0, so the coefficient at exp(/c0) 
in (|4.13p is negative. This contradicts the fact that I is a Laplace transform 
of a discrete probability measure (with non-negative point masses). Thus 
r € N. Concluding, we get (1). 

In order to get (2), let us now consider I with r = —1/A G 2N on = 
{9 : ai + a2exp(0) < 0} (condition l{9) > 0 is satisfied on 0 € D[). Since r 
is even, 

l{9) = (ai + a 2 exp(61))'’ = (—ai — 02 exp(6l))'’, 

and we analyze it on D[ = {9 : —ai — 02 exp(0) > 0}. Denoting di = —ai 
and 0.2 = — 02 -, we arrive at the case considered in (1). Therefore di, d 2 > 0, 
which yields a!i,a2 < 0. Thus (2) follows. □ 


Acknowledgement. I would like to thank Jacek Wesolowski for encour¬ 
agement and many helpful discussions. 


Appendix A. 

A.l. Proof of Proposition [421 In order to prove Proposition 14.21 we will 
analyze the roots of (12.31) and examine corresponding R functions. Our aim 
will be to eliminate the solutions of ()4.5I) that do not lead (via L = 
to Laplace transforms of probability measures. To do so, we will treat R as 
a function of 9i only (with fixed 02)) and we shall use auxiliary Lemma lA.ll 
and Lemma lA.21 presented below to reject some inadmissible solutions. 

Lemma A.l. Let r > 0. Assume that Pm is a polynomial of degree m over 
M and define 

k 

f{9) = Pm{d) + ^Ajexp(Aj0) + B9exp{'j9), 0 G 0, 

i=l 


(A.l) 
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where 0 contains some neighborhood of zero, Ai < A2 < • • • < Afc (and none 
of them is zero), and 7, B, Ai are some real numbers. 

If f^ is a Laplace transform of a probability measure then Pm = Pq and 
B = 0. 

Proof. Since / is defined in a neighborhood of zero, we consider a charac¬ 
teristic function <f>{t) = f^{it). We have 

k 2r 

Pm{it) + -Bit exp(7zt) -|- Ai exp(Ajit) 

i=\ 

Function \4>\ is bounded on M as the absolute value of a characteristic func¬ 
tion of a probability measure. Therefore Pm = -Bo and B = 0. □ 

Lemma A.2. Let r > 0 and 
1 

/(^) = X] cos(7j 6') -b {Aij + OBij) sin(7j6»)] 

j=o 

3 

-b ^ Aj exp(Aj0), 9 € Q, 
i=2 

where 0 contains some neighborhood of zero, be a function with all the pa¬ 
rameters being some real numbers. Furthermore, assume that Aq -b *70 7^ 
Ai -b i7i. Then if f^ is a Laplace transform of a probability measure, then 
Aoj = Aij = Bqj = Bij = 0, j = 0,1. 

Proof. Let 4>{t) = f^{it) be a characteristic function corresponding to 
then 



\m" = 

( e (^oj + itBoj) ---b e {Aij + itBij)- 


j=0 


2i 


+ A 2 exp(iA2t) -b A3 exp(iA3t) 


2r 


Since \(f \ is bounded on M, we arrive at the conclusion. □ 

Proof of Proposition 14.31 We shall separately consider all situations regard¬ 
ing the roots of ()2.3p . 


A. 1.1. Four single real roots of (12.31) . Let Ai, A2, A3, A4 be the roots of (j2.3p . 
Then R is of the form 

4 

R{9i,92) = y]] Ai((92)exp(Aj6»i), 
i=l 


(A.2) 
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where Ai{-), i = 0,...,3 are some real functions. Plugging (jA.2p into (j4.4p 
we obtain 

A2 a*( 02) - a\iAi{e2) - bA'M - eAAi{e2) = 0, i = 1,2,3,4. 
Therefore the explicit formulas for j4j(-)’s are 

Ai{92) = Ai exp ^ — —0^ , i = 1, 2, 3,4, 

and Ai, i = 1,2,3,4 are some real constants. Hence 

(A.3) R{ei,e 2 ) = Y.Ai exp . 

A.1.2. One double and two single real roots of (12.3h . Let Ai,A2 be single 
roots and A3 a double root of (12.3p . Then R is of the form 

2 

(A.4) R{9i,92) = ^ Ai{92) exp(Ai 6 'i) + (^ 3 ( 6 ( 2 ) + A 4 ( 6 ' 2 ) 6 'i) exp{Xs9i), 

2 = 1 

where Aj(-), i = 1, 2, 3,4, are some real functions. Since R~^^^ is a Laplace 
transform of a probability measure, using Lemma lA.ll we conclude that 
A4 = 0. Plugging (IA.4P into (14.3p we obtain 

AfAi(02) - a\A(9) - bA[i92) - eAAi{92) = 0, i = 1,2,3. 

These yield 

Ai{92) = Ai exp ^ — —92^^ , i = 1,2,3, 

and 

R{9i,92) = ^ Ajexp ^Aj^i + ^^ — —92^] ■ 

A.1.3. One single and one triple real roots of (12.31) . Let Ai be a single and 
A2 a triple real root of ()2.3p . Then 

-^(^1) ^2) = ^1(^2) exp (Ai^i) + (^2(^2) + ^3(^2)^! + ^4(^2)^!) exp(A20i), 

where Aj(-), i = 1, 2,3,4, are some real functions. Analogously as in the pre¬ 
vious subsection, using Lemma [A.II (for 6*2) exp(—A20i)), we conclude 
that A3 = A4 = 0. Function R in such case is of the form 

R{9i, 92) = '^Ai exp ^Ai^i + ^^— —92^] , 

where Ai, i = 1,2 are some real functions. 
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A. 1.4. Quadruple real root of (j2.3p . Let Ai be a quadruple real root of (j2.8p . 
Then 

R{ 0 i, 62 ) = Ai exp ^Ai 6 »i + — — ^ — — 02 ^ , 
where Ai is a real constant. 

A. 1.5. Two double real roots of (12.31) . Let Ai and A 2 be two distinct double 
roots of (12.3p . Steps analogous to the ones taken in Section lA.1.21 with the 
help of Lemma lA.ll yield 

R{ 0 i, 62 ) = '^Ai exp rAi 6 »i + ^ — — 02 ] , 

2=1 ^ ^ 

where Ai, i = 1 , 2 , are some real constants. 

A. 1.6. Two single real and two single complex roots of (12.31) . Let Ai, A 2 be 
two real and A 3 + ^ 73 , A 3 — 273 be two complex roots of (j2.3p . Then R takes 
the form 

R{0u92) = J^Ai(02)exp (xiOi + ^ 

i=l ^ ' 

exp(A 36 (i) (^ 3 ( 02 ) cos( 736 (i) + ^ 4 ( 6 ( 2 ) sin( 736 'i)). 
From Lemma lA.21 we conclude that A 3 = A 4 = 0, hence 

R{0i,02) = ^ Ajexp ^Aj^i + ^——02 

i=l ^ 

where Aj, i = 2, 3 are some real constants. 

A.1.7. Four distinct complex roots of (12.3p . Let Ai + ^ 71 , Ai — ^ 71 , A 2 + ^72 
and A 2 — *72 be the roots of (12.31) . Then R takes the form 

R{ 0 i, 02 ) = exp(Ai 6 'i) (Ai( 6 ' 2 ) 003 ( 716 ( 1 ) + A 2 { 02 ) sin( 7 i 6 (i)) + 

exp(A 26 (i) {A^{ 02 ) 003 ( 726 ( 1 ) + Ai{ 02 ) sin( 726 (i)). 

Using Lemma [A. 2 1 we conclude that there does not exist probability measure 
with the Laplace transform L = R~^l^ . 

A. 1.8. Two double complex roots of dOP. Analogously to the preceding 
case, there is no probability measure corresponding to R 



R{01,02) 


exp(Ai6(i) 


COs(7i6(i)(Ai(6'2) +0iA2{02)) + 


sin(7i6(i)(A3(6(2) + 0iAi{02)) , 


where Ai + *71 and Ai — *71 are double complex roots of ()2.3I) and Aj, 
i = 1, 2, 3,4, are some real functions. 
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A.1.9. Double real root and two single complex roots of (12.,ip . Let Ai, \ 2 +i ^2 
and A2 — i ^2 be the roots of (j2..Sp , then 

R{0i,02) = exp(Ai(9i) (Ai(( 92) + 0iA2{92)) 

+ exp(A 26 'i) (^ 3 ( 02 ) 003(726(1) + ^4(^2) sin(726(i)). 

From Lemma lA.ll and Lemma lA.21 we obtain A2 = A3 = A4 = 0. Hence R 
takes the form 

i?(0i,02) = exp (^Ai0i + . 

The conclusion of Proposition 14.21 follows by a straightforward aggregation 
of the above points. □ 
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